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Abstract. We study two classes of linear difference differential equations 
analogous to Euler-Cauchy ordinary differential equations, but in which mul- 
tiple arguments are shifted forward or backward by fixed amounts. Special 
cases of these equations have arisen in diverse branches of number theory and 
combinatorics. They are also of use in linear control theory. Here, we study 
these equations in a general setting. Building on previous work going back to 
de Bruijn, we show how adjoint equations arise naturally in the problem of 
uniqueness of solutions. Exploiting the adjoint relationship in a new way leads 
to a significant strengthening of previous uniqueness results. Specifically, we 
prove here fTheorem llll that the general Euler-Cauchy difference differential 
equation with advanced arguments has a unique solution (up to a multiplica- 
tive constant) in the class of functions bounded by an exponential function 
on the positive real line. For the closely related class of equations with re- 
tarded arguments, we focus on a corresponding class of solutions, locating and 
classifying the points of discontinuity. We also provide an explicit asymptotic 
expansion at infinity. 



1. Introduction 

We study two classes of linear difference differential equations analogous to Euler- 
Cauchy ordinary differential equations, but in which multiple arguments are shifted 
forward or backward by fixed amounts. Special cases of these equations have been 
studied quite extensively in the context of certain problems in number theory and 
combinatorics — see 32 for a brief overview. Here, we take a broader viewpoint. As 
in [S], we focus primarily on the advanced-argument linear difference differential 
equation 

m 

uq'{u) — ''^^ajq{u + Vj), u > 0, (1-1) 

3=0 

in which the argument u is incremented by the non-negative real numbers = vq < 
vi < ■ ■ ■ < Vm, and the coefficients ao, ai, . . . , am are arbitrary complex numbers. 
Here, however, we study (If in conjunction with the related delay differential 
equation 

rn 

{up{u)y = -^^ajpiu-Vj). (1.2) 
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It turns out that by judiciously pairing equations from each class, one can infer 
properties of solutions to one equation from properties of a solution to the other. 
We exploit this relationship to significantly strengthen previous uniqueness results 
concerning equation Our main uniqueness result is that Hl.l|l has a unique 

solution (up to a multiplicative constant) in the class of exponentially bounded 
functions. 

To exploit the relationship between the two equations 1)1.1(1 and ((1.2|l . we need 
to amass a certain amount of information concerning a solution to the latter. In 
this direction we establish new results, including the nature and location of discon- 
tinuities, and an explicit asymptotic expansion. 

The organization of the paper is as follows. In a brief historical overview is 
provided. We discuss uniqueness of solutions to ((1.1(1 in the class of polynomially 
bounded functions in ^ The next section focuses on a special case of ((1.1(1 arising in 
sieve theory and the study of cycle lengths of a random permutation. The particular 
properties of the parameters in this case are exploited to show that any solution 
to the underlying equation which is not polynomially bounded must be very wild 
indeed. The techniques of 21 prepare the ground for our exponential uniqueness 
theorem. We motivate the introduction of the adjoint equation in 50 Information 
concerning the asymptotic behavior of solutions to the adjoint equation is needed, 
and the relevant results are presented in ^ Details of the proof of our exponential 
uniqueness theorem are relegated to the penultimate section. 

Notation. The Bachmann-Landau 0-notation and the standard notion of asymp- 
totic equivalence are briefly recalled in the forms that we wish to use them here. If 
g is a positive function of a positive real variable, the symbol 0{g(u)), < u — > oo 
denotes an unspecified function / for which there exist positive real numbers uq 
and B such that |/(w)| < Bgiu) for all real u > uq. Two functions / and g are 
asymptotically equivalent at infinity, written f{u) giu), < u — s- oo, if for every 
e > 0, there exists a positive real number uq such that \f{u) — g{u)\ < e\g{u)\ 
whenever u > uq- 

2. Historical Overview 

Ever since Dickman [211 showed that the asymptotic density of the positive 
integers n < x having no prime factors exceeding x^/" is given by the continuous 
solution p{u) to the difference differential equation 

up'{u) = —p{u — 1), u> 1, 

with boundary condition 

p{u) = 1, < u < 1, 

it has been known that difference differential equations arise in the study of certain 
problems in analytic number theory. We now know of at least four types of number 
theoretical problems in which difference differential equations arise. These may 
be loosely classified according to the context of the original problem. The first 
type, "psixyology" (37j encompasses those problems associated with determining 
the probability distribution of the prime factors of an integer [HI IllL 1121 114L 1241 [301 • 
The second type, "sieve theory" [2H1 is concerned with the problem of estimating the 
size of a finite set of integers after certain residue classes have been eliminated [^ 
[ini[iniliniinilIHllinil23l2ZilSni- The third type deals with the problems relating 
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to the determination of the cycle decomposition of a random permutation |2 13 1251 
EEl 1^ OHl OH J and is closely related to the first P . The fourth type focuses on 
estimating incomplete sums of multiplicative functions pi ll3[l^ . 

Despite the extensive literature on the subject, apart from the treatment of a 
few isolated cases |3|H1ETIE2I2S1EI1E11ES1I1I1I111 the main focus has been 
on the application to problems in number theory and combinatorics as opposed 
to the study of the underlying difference differential equations as a problem in its 
own right. Although a systematic study of equation Hl.l|l in full generality was 
initiated in 9 , there is still additional work that remains to be done with regard 
to providing a comprehensive treatment. The present paper can be viewed as a 
further contribution in this respect, in that we explore more deeply the question 
of uniqueness, and reveal more fully the role of the adjoint equation (|1.2(l . It is 
hoped that some of the analytical techniques used herein — and more generally our 
ongoing program of study — will be of interest to researchers in number theory and 
the wider differential equations community as well. 



3. POLYNOMIALLY BOUNDED SOLUTIONS 

As in [HI, we fix a non- negative integer m, real numbers = vq < vi < ■ ■ ■ < Vm, 
complex numbers ag, ai, . . . , a™, and put (3 := Without loss of 

generality, we assume aj ^ for 1 < j < m. Define a function q* on the positive 
real half- line (0..cx)) as follows. 

Definition 1. For w > 0, let 

q*(u) := — — — / a;~''~^exp< —ux+y aj I dt>dx, (3.1a) 

n-P) Jo I fr{ Jo t 1 

if /3 is a negative integer, and 

, r(/3 + i) f _g_^ r ^ re"^*-ij. 

o iu) := z ' exp < uz — > a, / at 

otherwise. 

The contour in (|3.1b|) starts at — cx), hugs the lower side of the negative real axis, 
then circles the origin in the positive (counter-clockwise) direction before returning 
to —OQ along the upper side of the negative real axis. By deforming the contour, 
one can readily check that (|3.1a|) and Ij3.1b|) are equivalent if 5R(/3) < and /3 is not 
a negative integer. It follows [5| Cor. 2.5] that if aj is real for each j = 0, 1, . . . , to 
and /3 < 0, then q* is positive, convex, and log-convex. If /3 is a non-negative 
integer, then we can close the contour and deform it into the unit circle, so that 
then q* is a polynomial of degree (3. Otherwise, the singularity of the integrand is 
a branch point at z = 0, and q* is a transcendental function. A detailed treatment 
of the various properties of q* , including asymptotic expansions, behavior at 0, and 
a representation of q* [u) as an exponential of a Hellinger type integro-differential 
operator acting on the monomial is given in 9 . 

The main interest in q* stems from the fact that it provides a C°° solution to the 
difference differential equation (|1.1|) . More specifically, combining Proposition 2.2 
and Corollary 2.4 of yields the following result. 



\dz, (3.1b) 
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Proposition 1. The function q* of Definition^ satisfies the difference differential 
equation Hl.l|l and the asymptotic formula 

q*{u)^u^, 0<u^oo. (3.2) 

We recall jHl Theorem 3] that there is at most one function q satisfying the differ- 
ence differential equation (|f and the additional condition that q{u) is asymptotic 
to a fixed power of u for large real values of the argument u. Combining this result 
with Propositionn]yields the following uniqueness result, which we record for future 
use: 

Theorem 2. Suppose that q satisfies the difference differential equation Hl.l|l and 
the additional requirement that there exists a complex number t such that 

q{u) ^ u'^ , < M — > oo. 

Then t = [3 and q = q* , where q* is as in Definition^ 

In other words, the function q* of Definitionn]is the unique solution to the differ- 
ence differential equation 1)1. l|l in the class of functions asymptotic to a monomial 
at infinity. Theorem 4 of strengthens Theorem|21by widening the class of admis- 
sible functions to those functions which are polynomially bounded (i.e. majorized 
by a polynomial) in the right half-plane. The proof in (2j of this latter result is un- 
necessarily complicated, applying residue theory to the inverse Laplace transform of 
a suitably high order derivative of a supposed solution. Before proceeding further, 
it may be of interest to give a simpler proof which has the additional advantage 
of requiring only that the supposed solution be majorized by a polynomial on the 
positive real half-line. 

Theorem 3. Suppose thatq satisfies the difference differential equation Sup- 
pose further that there exists a real number r such that 

q{u) ^ 0(u''), < M ^ oo. (3.3) 

Then there exists a complex number A such that for all u > 0, q{u) — Aq*(u), where 
q* is given by Definition^ The multiplicative constant A may of course depend on 
the parameters Uj, Vj {j = 0, 1, 2, ... , m) but not on u. 

Proof. For any complex number A, Theorem |21 implies that q — Aq* is the unique 
function satisfying both the difference differential equation and the asymptotic 
formula 

q{u) - Auf^, < u ^ oo. (3.4) 

Accordingly, it is enough to show that our hypotheses imply the asymptotic for- 
mula H3.4)l holds for some complex number A, or equivalently, that 

lim u^^q{u) 

u — *oo 

exists. 

We can write the difference differential equation in the form 

m 

(u^'^q{u))' = u"^"^ ^ Oij{q{u + Vj) - q{u)}, u> 0. 
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We have dropped the term in the sum that corresponds to j = 0. It vanishes 
because vq = 0. It foUows that for u > 0, 

/It m 
t-^-^ aj{q{t + Vj) - qit)}dt 

= g(l) + Vaj / t^'^^M q'{x + t)dxdt. (3.5) 

j=l "'1 "'0 

Let c := 3?(/3), and suppose that the polynomial growth requirement H3.3|l holds. 
Since increasing r only weakens the hypothesis, there is no harm in assuming that 
r = n + c+1/2 for some non-negative integer n. The stipulation that r — c lies 
halfway between two consecutive integers is not strictly necessary; however, as it 
precludes the possibility that an integer power of m multiplied by u^~'^ will integrate 
to a logarithm as opposed to a power of u, the number of cases to be considered in 
our subsequent estimates is conveniently reduced. 

By the difference differential equation and our hypothesis, we have 

rn 

\uq'{u)\ < |aj||(7(it + Vj)\ = 0{u^), < u ^ oo, 

and hence q'{u) — 0{u^^^). Using this bound in the integral 1)3.5(1. we conclude 
that 

u-f^qiu) = q{i) + o r^^h"--^ d?j = 0(1) + o(m"-"-1), 

and hence q{u) = 0{u'^) + 0{u^^^). By repeating the previous steps if necessary, 
we arrive at the growth estimate 

q{u) = O(m^), < u — > oo. 

Using this latter bound in the difference differential equation Hl.lfl yields 

m 

\uq'{u)\ < ^ |aj||q(M + Vj)\ — O(w^), < m ^ oo, 

3=0 

and hence q'{u) — 0{u'^^^). It follows that the integrals 

/ q\x + t)dxdt (j = 0, 1,2, . . . ,m) 



all converge. Therefore, we may let u — > +oo in H3.5|l . from which we infer that the 
limit 



^ /'OO nVj 

A:= lim u-'3q(w) = g(l) +y a, / rl^'^ / q' 



{x + t) dx dt 



exists. In other words, 

u — > oo, 

and by our initial remarks, the proof is complete. □ 

In ij7| we show that the conclusion of Theorem 13 holds with H3.3|l replaced by 
the weaker hypothesis that q be majorized by some fixed power of an exponential 
function. Before proving this stronger result in full generality, it is instructive to 
examine an important special case which is used in the proof of the general case. 
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4. The Equation — Kq{u) - nq{u + 1) 

Let K be a positive real number. The difference differential equation 

{uq{u))' ^ Kq{u) - Kq{u + 1), u > 0, (4.1) 

was introduced by Iwaniec [33] as an adjoint equation for the upper and lower 
bounding sieve functions of Rosser's sieve, and is a special case of with to = 1, 
vi — 1, a{) — K — 1, ai = —K and (3 = + ai = —1. The role of k in sieve 
theory is to measure the average number of residue classes being deleted for each 
prime used in the sifting. Equation (|4.1|l and its close cousin with the minus sign 
replaced by a plus sign have been studied fairly extensively in the context of sieve 
theory |3|Hl[n|[THl[T3|iniEIlEllEll- The case k = 1 of glj also occurs in the 
problem of determining the asymptotic average size of the largest prime factor of 
a random integer and the longest cycle of a random permutation |25l 1261 1841 1^ . 
More specifically, the average cycle length of the longest cycle in a permutation on 
n symbols |40| and the average number of digits in the largest prime factor of an 
rt-digit number |34| are both asymptotic to nq[{l) / qi{l) , where 

qi{u) J exp |— ua; — J dt j dx, u > 0, 

satisfies l|^?T|) with k = 1. We also have gllllll 

#{n eZ:l<n< x,P2{n) < Pi(n)i/"} = xe'<qi{u) + 



. log X ^ 

where Pj{n) denotes the jth largest prime factor of n and 7 is Euler's constant. 

By the results of ^ we know that the polynomially bounded solutions to (|4.1() 
are all constant multiples of q*, where here 

q* (u) — J exp ■j^ — ua; — K J dtjdx, u>0. (4-2) 

It turns out that any solution of (|4.1() which is not a constant multiple of H4.2|l must 
oscillate wildly. More precisely, we show that if q satisfies l)4.1fl with k > 0, and if 
there is no constant A for which q = Aq*^ then q(u) ^ 0(exp(Aw)) for any fixed 
A > 0, and for all u > 0, g changes sign infinitely often in [u, 00). To prove this, we 
first establish the following special case of our main uniqueness result. 

Lemma 4. Suppose that q satisfies the difference differential equation with 
/3 = — 1. Suppose further that there exists a positive real number A such that 

q{u) = 0(e^"), < u 00. 

Then there exists a complex number A such that for all u > 0, q{u) — Aq*{u), 
where q* is given by Definition^ with f3 = —1, i.e. 

f°° ( 1 — e~'"j* 1 
q{u) ~ ^ J — ux + ctj J dt j dx, u > 0. 

Proof. Partial summation enables us to rewrite the differential difference equa- 
tion in the form 

m 

{uq{u)y ^ Cj{q{u + Vj) - q{u + Vj-i)}, u > 0, 
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where 

Now integrate. There exists a complex number A such that for all u > 0, 

uq{u) = A + ^Cj q{u + 1) dt. (4.3) 



Let 

A 



The exponential growth requirement implies that there exist uq > M and B > Q 
such that for all u > uo, \q{u)\ < Be^^. Inserting this inequality into 1)4. 3() yields 



\uq{u)\ < \A\+bY,\cj\ / 

.7 = 1 -^"^-1 



;^("+*) dt, u > Mo- 



It follows that 

r)\ < . 

u Xu 



g(«)i< V + ^Ei^^> ^ =^ V + (^-^^ 

Now we claim that \q{x)\ < 2\A\/x for x > uq. To prove the claim, fix a; > mq and 
observe that if i?e'^'^/4 < \A\/x, then 14.4|l gives \q{x)\ < 2\A\/x. On the other hand, 
if Be^'^/A > \A\/x, then we must have i?e^"/4 > for all u > a;. Since (|4.4|) 

holds for all u > uq, we then get 

< , u>x. 

We can now insert this latter inequality back into (|4.3|l and repeat the previous 
reasoning with B replaced by B/2. In general, if n is a positive integer such that 

Be^^ 1^1 
4 • 2"-i ^ 

then iterating the previous argument n times will yield 

, , U| Be^" , , 

|'ZHI<— + ^^, u>x. 4.5 
u 4-2" 

Note that the bound ()4.5|l is valid for u > a;, where x is the same as above. If 
A 7^ 0, let n be the least positive integer such that 

Be^^ 1^1 
4 ■ 2" ~ 

Then (|4.5|l gives 

\q{x)\ < (4.6) 

Since x > uq is arbitrary, (|4.t)|l must hold for all x > uq. On the other hand, 
a A — 0, then letting n — > cx) in H4.5|l shows that q{u) — for all u > x. But 
again, since a; > uo is arbitrary, we must have q(u) = for all u > uo- In other 
words, (|4.6|l also holds for all x > uq when A = 0. This establishes the claim. 
Invoking Theorem |31 now completes the proof. □ 
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We can now show that any solution to H4.1|l with k > must be wildly oscillatory. 

Theorem 5. Suppose that q is a solution to the difference differential equation (|4.1() 
with K > 0. Suppose further that q is not a constant multiple of q* , where q*{u) is 
given by (|4.2I) . Then q{u) ^ 0(exp(Au)) for any fixed A > 0, and for all u > 0, q 
changes sign infinitely often in [u,oo). 

Proof. Note that /3 = — 1 in H4.1|l . The conclusion q{u) ^ 0(exp(A-u)) for all 
A > is an immediate consequence of Lemma 0] In particular there exists a 
strictly increasing sequence of positive real numbers ui < U2 < . . . satisfying 
lim„^oo M„ = oo and |g(M„)| > exp(u„) for all positive integers n. Integrating the 
difference differential equation H4.1|) yields 

i-u+l 

uq{u) = A~ K q{t) dt, u > 0, (4.7) 



where A is a constant of integration which may depend on k, but not on u. There 
is a positive integer N such that if ti > iV then u„ > 1 and exp(M„) > \A\. Let 
n > N and set u = u„. Suppose that q{u) > 0. Then (|4.7|) yields 

ti+i 

q{t) dt^ A- uq{u) < \A\ - e" < 0. 
Since k > 0, it follows that 

q{t) dt < 0, 

and hence q must change sign from positive to negative somewhere in the half-open, 
half-closed interval {u, If, on the other hand, q{u) < 0, then H4.7|l yields 

u+1 

q{t) dt = A- uq{u) > -\A\ + e" > 0, 

and hence q must change sign from negative to positive somewhere in (m, u -|- 1] . 
Thus, we have shown that q has a sign change beyond m„ for each n > N, and the 
proof is complete. □ 

5. The Adjoint Relation 

We'd like to remove the restriction f3 — —1 in Lemma 01 A careful examination 
of the proof of Lemma ^ suggests that the iterative argument employed therein 
succeeded because the lengths of the ranges of integration in the representation (|4.3() 
are all independent of u. Note that H4.3|l is equivalent to 



m .u 

uq{u) = A + aj I 

1 = 1 Ju-Vj 



q(t + Vj) dt. 



In order to recover this property in the case when /? — 1 , let us seek a function p 
such that for some constant A, the equation 

m „„ 

up{u)q{u) = A + '^aj p{t)q{t + vj) dt (5.1) 

, = 1 Ju-Vj 
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holds for all sufficiently large values of u. Let us suppose further that beyond some 
point, p is differentiable. Differentiating (|5.HI with respect to u reveals that 

m rn 

uq'{u)p{u) + up'{u)q{u) + p{u)q{u) — p{u) ^ ajq{u + vj) - q{u) ^ ajp{u — vj). 

If we now substitute the right hand side of the difference differential equation 

for uq'{u) in this latter equation, the ffi'st sum on the right hand side can be 

cancelled with the corresponding sum on the left. The result is 

m 

aop{u)q{u) + up{u)q{u) + p{u)q{u) = -q{u) ^ ajp{u -vj), 

which is certainly the case if p satisfies the delay differential equation (|1.2I) . Con- 
versely, it is clear that if q satisfies Hl.lll and p satisfies (|1.2f) . then (|5.1I) must also 
hold for all sufficiently large values of u. 

6. The Function u i-^ p{u, a, b) 

Our plan is to use the adjoint equation (|5.1|l to deduce the asymptotic behavior 
of q{u) for large u using only a very weak estimate on the growth rate of q{u). 
In order to do this, we need reasonably precise knowledge concerning the rate of 
growth of a non-trivial solution p{u) to the delay differential equation 1)1. 2|l as u 
increases without bound. Most of what we need can be found in |35) . where the 
existence of an asymptotic expansion for p for the case (in our notation) K(ao) < 
is proved, and the first term of the expansion is determined explicitly. Nevertheless, 
it is worthwhile to obtain some additional results, and for this it is useful to extend 
Wheeler's ^202] more comprehensive treatment of the m — 2 case to arbitrary m 
and complex aj. Therefore, with regard to the choice of boundary condition and 
the location and classification of discontinuities, our development in this section 
more closely parallels that of j42| . However, to provide an explicit formula for all 
coefficients in the asymptotic expansion, we find it more convenient to extend the 
technique used in . To facilitate comparison with our results, we set a = 1 -|- ag 
and let b denote the vector (ai, a2, . . . , am)- 

Let 

m 

where, as usual. 



7 = lim 



.7 = 1 ■' ^ 



is Euler's constant. Consider a particular solution p(u) = p{u,a,b) of the delay 
differential equation H1.2I) . which is defined uniquely for all real u by the following 
six conditions: 

m 

{up{u)y = -^ajp{u-Vj); (6.1a) 

p{u) — whenever u < 0; (6.1b) 

if 0<^<i;i; (6.1c) 

r(l - a) 
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p{u) is continuous for u > when 5R(a) < 1; (6. Id) 

p{u,a+ l,b) — -^p{u,a,b); (6.1e) 
au 

p is continuous from the left at every point. (6. If) 

We need to prove that a unique such function exists. First assume 5ft(a) < 1. The 
function defined by (|6.1c|l is integrable for < u < vi and the delay differential 
equation Ht).la|) can be rewritten as 

m 

(uXu))' = -u"--^ ajp{u - Vj), 
i=i 

which can be integrated forward on successive intervals. That is, for each positive 
integer n, ifn<u<n + l then 



p{uvi) 



-""In^pinvi)- I t''-^Yajp{tvi-Vj)dt\. (6.2) 



This establishes uniqueness and every condition except H6.1ef) . To see that (|6.1ef) 
holds, simply differentiate (|6.1a() and (|6.1c|l and apply the functional equation for 
the gamma function. For the case 5R(a) > 1, we use the fact that if n is a positive 
integer, then 

d 



p{u,a,b) = p{u,a-n,b), 

which follows from (|6.1e|l . 

From H6.2|l , it follows that p has a two-sided derivative at every point x for which 
p is continuous at each of the points x — Vj, 1 < j < m. Otherwise, we need to take 
derivatives from the left in H6.1al) and (|6.1e|l . This we can do, again by 1)6.2(1 . since 
p is continuous from the left at every point. 

Theorem 6. 

(i) // 5R(a) is not an integer, then the set of points at which p is discontinuous 
is precisely {nvj : n G Z, < n < 5R(a), 1 < j < m}. At these points, p has 
a finite limit from the left, and is unbounded from the right. Moreover, if n is a 
non-negative integer, 3ff(a) > n, and 1 < J < n-, then 

^(")~n!.^"r(n-a+l)("-""^) ' ^ ^ ^-^+' 
whereas lim p{u) ~ p(nvj) is finite. 

u — >nvj — 

(ii) // 5ft(a) is a non-negative integer, but a is not itself an integer, then the set 
of points at which p is discontinuous is precisely {nVj : n G Zi, < n < 5R(a), 1 < 
J < m}. The discontinuities consist of two types. If n is a non-negative integer, 
5ft(a) > n, and I < j < ni, then as in (i), 



whereas lim p{u) — p{nvj) is finite. If n — 5R(a), then p is bounded in a neigh- 

u — >nvj — 

bourhood ofnvj, but lim p{u) does not exist. 

u — >nvj + 

(iii) If a is a positive integer, then p has only finite jump discontinuities. These 
can only occur at the points nVj, where n is an integer, 1 < n < a, and 1 < j < m. 
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(iv) Finally, p{u, 0, b) has a finite jump discontinuity at u = and no other 
discontinuities. 

Sketch of Proof. The statements of the Theorem are true if 5ft(a) < 1 as we aheady 
noted. The other cases can be estabhshed by induction, using the equation 

d a ^ m 

p(u, a+ l,b) — —p(u, a, b) = a, 6) > ajp(u — Vj, a, 6), (6-3) 

du u u 

which is certainly vahd if u is not an integer multiple oivj, 1 < j < m. For details, 
see 132 Theorem 1] □ 

Remark. The situation in (iii) is somewhat unsatisfactory, as in general all we can 
say is that if a is a positive integer and u is a point of discontinuity, then necessarily 
u — nvj for some positive integer n < a and 1 < j < m. Although it is typically 
the case that all such points are indeed points of discontinuity — for example, this 
is indeed true if the m positive real numbers vi, . . . , Vm are linearly independent 
over the integers — it may happen that for certain values of the parameters aj , Vj , 
the jumps arising from different terms in 1)6.3(1 cancel each other. As an example, 
consider the case m = 2, a = 2, 6 = (ai,Q;2) — (1,-2), vi = I, V2 = 2. Then 
Co = Ae"'. We find that the jump at 2vi+ is exactly cancelled by the jump at V2 + , 
so that p(2+, 2, b) = e'' = p(2— , 2, b). As noted by Wheeler [321, such cancellations 
cannot occur when m = I. 

We now turn to the problem of determining the behavior ofp{u) for large positive 
real u. We take the following result of Levin and Fainleib [SS] as our point of 
departure, recasting it in our notation, and with our boundary condition (|6.1c|l etc. 

Lemma 7 (jHSli Lemma 1.3.1). Let di{a) < 1, and let (p{u) = {ulogu)/v,n + 0{u), 
< u — > oo. We have the following behavior for p{u) as < u — s- cx). 

(i) If [3 is a non-negative integer, then p{u) — 0(exp(— (/?(!*))). 

(ii) If (3 is a negative integer, there exists a polynomial rp of degree —[3 — 1 such 
that p{u) = rp{u) + 0(exp(— 

(iii) For any complex (3 not an integer, there is an asymptotic expansion p{u) ^ 

Levin and Fainleib '35' did not identify the polynomial rp in (ii). In addition, 
they gave only the first term in the asymptotic expansion (iii). We shall determine 
all the coefhcients c„ and the polynomial rp explicitly. We shall also remove the 
restriction on 5R(a). To carry this out, we need the following formula for the Laplace 
transform of p. 

Theorem 8. Let 5R(a) < 1 and 5R(s) > 0. Then 

J e"^Xu) d-u = s'^ exp I - 2J y dtj. 

If P is a non-negative integer, the formula is valid for all complex numbers s. 

Theorem|Slis actually just a restatement in our notation of a corresponding result 
derived in |SS|- Although the region of convergence was not discussed there, both 
the existence of the Laplace transform and its region of convergence follow easily 
from the growth behavior of p as given in Lemma H 
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Let n be a non-negative integer. Recall the polynomial equation (3.2)] 

d_ 

dz 



e^J* - 1 

Qn{u,b) 

z=0 



Remembering that b denotes the vector of coefBcients (ai, . . . , am), we let Qn(u, —b) 
be the polynomial obtained by negating each aj {I < j < m) in Qn{u, b). With this 
notation, we can now restate Lemma [3 with the coefBcients c„ and the polynomial 
explicitly identified. 

Theorem 9. Let (p{u) be as in the statement of Lemma^ As < u — > oo, we 
have the following behavior for p{u) . 

(i) If (3 is a non-negative integer, then p{u) = 0(exp(— (y5(w))). 

(ii) If P is a negative integer, then 

Pin) . ^ .yrA, --^-^-" + 

n=0 • V / J- 

(iii) For any complex number j3 not an integer, we have the asymptotic expansion 

y- (-i)"Qn(o.-^) ,-.^i^» 

2^ , r(-/3-n) 

Proof. In view of Lemma there is nothing to prove for (i), and the error term 
need not concern us in (ii). Fix a positive integer h. It is sufRcient to consider 
the case 5R(a) < 1 and K(/3) < \ — h, since as noted in [ESJ, the remaining cases 
reduce to the problem of differentiating asymptotic expansions, and by H6.1a|l this 
is legitimate. Let < e < 1 be such that — 5R(/3) — h — e > —1. Substituting 

h-l 

p{u) = CnU-^-^-" + Oiu-^^^^^-''-'), < M -> OO 

n=0 

into the Laplace transform for p, we find that 

/•oo h — 1 

/ e-™p(M)du = ^c„r(-/3-n)s''+" + 0(s^('3)+''+^-i), s^0+. 

•^0 „=o 

On the other hand, Theorem |S1 and the definition of the polynomials Qn(0,—b) 
imply that 

Comparing the two expressions for the Laplace transform, we infer that 

(-1)" Q„(0,-6) 
n! r(— /3 — n) ' 

which proves (iii). If /3 is an integer and f3 > —n, we see that c„ — 0. It follows 
that if /3 is a negative integer, then 

^-^ n\ I — n — n) 
as stated in (ii). □ 
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7. Exponentially Bounded Solutions 

In this section we prove that, up to a multiplicative constant, the difference 
differential equation (|l.l|l has a unique solution in the class of functions majorized 
by a function of exponential growth. We first establish an inequality for certain 
exponential integrals that arise in the proof. 

Lemma 10. Let u, r and A be real numbers satisfying A > r+ = max(r, 0) and 
u> 1. Then 

H"' dt < 



X-r+' 

Proof. Denote the integral by /. An easy integration by parts shows that 

I < ^e^^u-"- + ^ f e^H-^^^dt. (7.1) 



If r < 0, then the coefficient of the integral on the right hand side of ()7.1|l is 
negative, and so the corresponding term can be dropped, yielding 

1 

A^ 

Since r < implies ?-+ = 0, the claim is established in this case. On the other 
hand, if r > 0, then t^'' > t^^^^ for 1 < t < u and hence 

/ < ie^"u-'- + ^ f e^H-'' dt = ie-^"?/-"^ + ^I. 
A A ./i A A 



Since A > r > 0, it follows that 



A-r 

Since r > implies r+ = r, the claim is established in this case also. □ 

Theorem 11. Suppose that q satisfies the difference differential equation ()l.l|l . 
Suppose further that there exists a positive real number A such that 

(7(w) = 0(e^"), 0<u^oo. 

Then there exists a complex number A such that for all u > 0, q{u) = Aq*{u), where 
q* is given by Definition^ The multiplicative constant A may of course depend on 
the parameters aj , Vj {j = 0, 1, 2, ... , m) but not on u. 

Proof. First, suppose that /3 = n is a non-negative integer. The difference 
differential equation (|l.l|l implies that q'{u) satisfies (|l.lll with ao replaced by 
ctQ — l, and that q'{u) = O (^"^e'''"), < m ^ oo. Hence, we deduce that the n + l^* 
derivative of q satisfies the conditions of Lemma 0] Accordingly, by Lemma ^ and 
Proposition ^ there exists a complex constant A such that 

q^"+'^\u) Au^^, 0<M^oo. 

It follows that 

q'-'^^u) ^ Alogu, < u ^ CXO. 

Also, q^"-' satisfies the difference differential equation Hl.l|l with /3 = 0. Theorem^ 
now implies that g^"-* is a constant multiple of q*, where q* is given by Definition^ 
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with f5 = 0. It follows that A = and (7^"^ is constant. Hence, there exist complex 
numbers ao, ai, . . . , a„ such that for all u > 0, 

n 

More specifically, Theorem |2| implies that q — a„Q„, where 

d 



exp < uz — y Oij dt 



is the polynomial obtained by setting /3 = n in (j3.1b|) of Definition^and deforming 
the contour into the unit circle. See (2] for additional properties of the polynomial 

Qn- 

For the general case, we assume that /3 is a complex number, but not a non- 
negative integer. Let c = 3?(/3). Since increasing A only weakens the hypothesis, we 
may assume that K := X — max(c + 1, 0) > 0. Let 

1 R JUL 

By hypothesis, there exist ui > M and B > such that for all u > ui, \q{u)\ < 
Be^". Let p{u) = T{—(3)p{u, a, b), where a, b) is as in ^ Since (5o(0, —b) = 1, 
Theorem |51 implies that p{u) ^ u~^~^ as < m — > 00. Therefore, there exists a 
positive real number U2 such that for all u > U2, ^m""^"^ < < 2m^'^^^. There 

must also exist a complex number A and a positive real number U3 such that the 
adjoint relation (cf. I5.1|l 

up{u)q{u) ^ A + ^a-i p{t)q{t + v-i) dt (7.2) 

j = l Ju-Vj 

holds for all u > U3. Let M4 :— max(ui, U2 + Vm, ^3, 1 + Vm)- Then for all u > U4, 
by Lemma 11111 we have 

\up{u)q{u)\ < \A\ + 2B"^\aj\ / t~^~^e^^'+^^'> dt 

j = l Ju-Vj 
m 

< |A| + 2Bif-iu-"-i^|ajle-^("+"^) 

< \A\ + !^Bu-^-ie^". 

8 

Let A' = 2A. It follows that 

|g(M)| < + — Be^" < lA'w'^l + iBe^", u > u^. (7.3) 

iu * 

Now let Mo ■— max(u4,c/A), and fix a; > uq. If jBe^^ < \A'x^\, then H7.3|) gives 
\q{x)\ < 2\A'xP\. 

On the other hand, if jBe^"^ > jA'x'^l, then we claim that jBe^" > \A'uf^\ for 
all u > x. To see the claim, let f{u) :— log i? — log 4 + Au - log | A' | - 5R(/3) logw. We 
have f{x) > and f'{u) = A — c/u > for all u > x. Therefore, / is positive and 
strictly increasing on the interval [a;, 00), and this proves the claim. Thus, in this 
case, H7.3|l gives \q{u)\ < \Be^" for all u> x. Now insert this latter inequality back 
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into the adjoint relation (|7.2|l . i.e. repeat the previous argument with B replaced 
by B/2. We get 

\q{u)\ < \A'u^\ + ifie^", u>x. (7.4) 
In general, if n is a positive integer such that 

then n iterations of the preceding argument will yield 

\q{u)\<\A'u^\ + ^, u>x. (7.5) 
Therefore, if A' ^ 0, let n be the least positive integer such that 



< \A'x' 



./3| 



4-2^ 

Then (|7.5|l implies that 

\q{x)\<2\A'x^\. (7.6) 
Since x > uq is arbitrary, the bound (|7.6I) must hold for all such x. On the other 
hand, if A' = 0, then letting n ^ cxd in H7.5|l shows that q{u) — for all u > x. 
But again, x > uq is arbitrary, so we must have q{u) — for all u> uq. In other 
words, the bound H7.6|l holds for all a; > uo and all complex numbers A' . Invoking 
Theorem O completes the proof. □ 

8. Final Remarks 

The introduction of an adjoint equation as an aid to studying solutions of differ- 
ence differential equations goes back at least to de Bruijn's study of the Buchstab 
function 10 . There and subsequently (see eg. (361 1151 41, 42 ), the adjoint equa- 
tion is used to deduce information about a solution to a special case of (|1.2|l from 
a solution to the corresponding special case of (|l.l(l . In this paper, we turned the 
process around, deducing information about a solution to (jl.l(l from the adjoint 
equation and a solution to H1.2|l . 

It is interesting to compare and contrast the behavior of p and q near zero and 
infinity. As in SQlet p{u) — T{—(3)p{u,a,b), where p{u,a,b) is as in |S1 Assume 
that P is not a non- negative integer. Then Theorem |51 implies that 

p{u) ^ £(-1)" (^'^^^ ^) g„(0, -6)u-'^-i-", < w ^ oo, 

n=0 ^ 

whereas Theorem 5 of |H| states that 

If in addition, di{aa) < 0, then Theorem 7(ii) of [H] can be restated as 

whereas ()6.1c|l implies that 

p{u) = u-^°-'^P^ f[ (vje-^y^ Q<u<v,. 
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In particular, lim„^o+ up{u)q{u) = linii,^+oo up{u)q{u) = 1. 
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